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Transitionless-based shortcuts for the fast and robust generation of W states 
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We propose a scheme to generate W states based on transitionless-based shortcuts technique in 
cavity quantum electrodynamics (QED) system. In light of quantum Zeno dynamics, we first ef¬ 
fectively design a system whose effective Hamiltonian is equivalent to the counter-diabatic driving 
Hamiltonian constructed by transitionless quantum driving, then, realize the W states’ generation 
within this framework. For the sake of clearness, we describe two stale schemes for W states’ gener¬ 
ation via traditional methods: the adiabatic dark-state evolution and the quantum Zeno dynamics. 
The comparison among these three schemes shows the shortcut scheme is closely related to the other 
two but better than them. That is, numerical investigation demonstrates that the shortcut scheme 
is faster than the adiabatic one, and more robust against operational imperfection than the Zeno 
one. What is more, the present scheme is also robust against decoherence caused by spontaneous 
emission and photon loss. 
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I. INTRODUCTION 

Quantum entanglement is an intriguing property of composite systems. The generation of entangled states for two or 
more particles is not only fundamental for demonstrating quantum nonlocality [l, 2}, but also useful in quantum information 
processing (QIP) H9- For three-qubit entanglement, there are two main kinds of entangled states , the W states fi| and the 
Greenberger-Horne-Zeilinger (GHZ) states [2j. These two kinds of entangled states cannot be converted to each other by local 
operations and classical communications. In recent years, the W states attract more attentions because of its robustness against 
qubit loss and advantages in quantum teleportation ra. So far, lots of theoretical schemes have been proposed to generate 
W states in different systems via different techniques p7HT4l] . There are two techniques famous for their robustness against 
decoherence in proper conditions and have been widely used in QIP: one is named stimulated Raman scattering involving 
adiabatic passage (STIRAP) including their variants fl5l - fl9| . and the other one is Quantum Zeno dynamics (QZD) [20l - [24| . 
Generally speaking, the adiabatic passage technique is robust against variations in the experimental parameters and atomic 
spontaneous emission. To restrain the influence of photon leakage on the fidelity, a widely used way is choosing parameters to 
reduce populations of the intermediate excited states. However, such operation inevitably increase the operation time. As is 
known to all, using adiabatic technique (we name it “adiabatic scheme” for short in this paper) asks for an adiabatic condition 
that the change of a system’s Hamiltonian in time is managed to be slow to make sure each of the eigenstates of the system 
evolves along itself. Using QZD method, by contrast, might be faster than using adiabatic passage. But that depends, especially 
in multiparticle systems. Usually, in a scheme based on QZD (we name it “Zeno scheme” for short in this paper), we consider 
the system’s Hamiltonian as H = H 0 b s + KH me a. s , where H a b s is the Hamiltonian of the quantum system investigated, K 
is the coupling constant, and ff meas is viewed as an additional interaction Hamiltonian performing the measurement. When 
K —> oo, the system’s effective Hamiltonian is approximated as Hz = ^2 (K£ n P n + P n H 0 b s Pn), where P n is the nth eigenvalue 
projection of ifmeus with eigenvalue Similar to the adiabatic passage, there is also a limited condition in a Zeno scheme 
that limits the system’s speed: the Zeno condition K —» oo. It has been confirmed by lots of schemes that using QZD for QIP 
is usually robust against photon leakage but sensitive to the atomic spontaneous emission. Therefore, in order to restrain the 
influence of atomic spontaneous emission, some researchers introduced detunings between the atomic transitions to decrease the 
population of atomic excited states. That also inevitably increases the operation time. In addition, the operation time required 
in a scheme via QZD always needs to be controlled accurately, which increases the difficulty to realize the scheme in experiment. 
As we know, the operation time for a method is the shorter the better, otherwise, the method may be useless because the 
dissipation caused by decoherence, noise, and losses on the target state increases with the increasing of the interaction time. 
Many experiments also desire fast and robust theoretical methods because high repetition rates contribute to the achievement 
of better signal-to-noise ratios and better accuracy. 

Therefore, fast and noise-resistant generation of entangled states becomes a research hotspot in recent years, especially, after 
the technique named “Shortcuts to adiabatic passage” (STAP) p5l . [26j was proposed. This technique is related to adiabatic 
passage but successfully breaks the limit of the adiabatic condition. It describes a fast adiabatic-like process which is not 
really adiabatic but leading to the same final populations with adiabatic process. Newly, STAP has shown its charm in 
theory and experiment I27H55I1 . In 2014, by using transitionless tracking algorithm under large detuning condition , Lu et al. 
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FIG. 1: (a) The experimental set-up diagram, (b) The atomic level configuration for traditional schemes (STIRAP and QZD). 
(c) The atomic level configuration for APF system. 


proposed an effective scheme to implement fast populations transfer and fast maximum entanglement preparation between two 
atoms in a cavity [30j. The idea inspires that using some traditional methods to approximate a complicated Hamiltonian into 
an effective and simple one first, then constructing shortcuts for the effective Hamiltonian might be a promising method to 
speed up a system. Soon after that, Chen et al. first combined invariant-based inverse engineering with Zeno subspaces to 
construct shortcuts to perform fast and noise-resistant populations transfer for multiparticle systems [2Sj. In their method, they 
demonstrated that besides constructing STAP, slightly broking down the Zeno condition under certain conditions is another 
simple way to speed up the evolution. Soon after that, similar ideas with slightly breaking the Zeno condition down are rapidly 
used to perform fast and noise-resistant QIP [H, |33|, IHfil - tssl ] . 

Motivated by the above analysis, we discuss how to construct STAP to rapidly generate W states by using the approach 
of “transitionless tracking algorithm” in cavity QED systems. Different from ref. [32| which proposed a method through 
combining Lewis-Riesenfeld theory and Zeno subspaces to generate a JV-atom W state by N + 1 atoms, we do not need to 
abandon any atoms. An 7V-atom W is fast generated directly by N atoms in one step. In order to explain the charm by 
using STAP to generate W states, we first give a brief description about generation of W states via two traditional methods 
(STIRAP and QZD). Then, we propose the scheme by using transitionless tracking algorithm in detail. The comparison among 
these three schemes demonstrates that the shortcut scheme is not only faster than the adiabatic one, and more robust against 
operational imperfection than the Zeno one. What is more, this method might be promising when it comes to generation of 
multi-level and multi-qubit entangled states, i.e., the singlet states. 

The rest of the paper is organized as follows. In section II, we give the model of atom-cavity system and introduce two 
schemes to generate W states via traditional adiabatic method and Zeno method. Then in section III, we use the transitionless- 
based shortcuts method to propose a fast and noise-resistant scheme to generate W states. The conclusion is derived in section 
IV. 


II. THEORETICAL GENERATION OF W STATES IN A THREE-ATOM SYSTEM 


For simplicity, we assume that three A-type atoms (cu, < 12 , 03 ) are trapped in a cavity (c) as shown in Figs. |T] (a) and (b), 
each atom has an excited state |e) and two ground states |/} and | g). Atomic transition |/} |e) is resonantly driven by 

classical held S2(f), and the transition | g) |e) is coupled resonantly to the cavity with coupling A. Under the rotating-wave 
approximation (RWA), the interaction Hamiltonian for this system reads 
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H I = Y, n k (t)\e) k (f\ + \ k a\e)k(g\ + H.c., (1) 

k=i 

where a is the annihilation operator of the cavity. We assume the initial state of the system is \ f, g, < 7 ) 1 , 2 , 3 10) c , the system will 
evolve within a single-excitation subspace V spanned by: 

l^i) = l/,5>fl)l,2,3|0) c , 

1^2) = |e,P,s)i,2,3|0) c , 

1^3) = |fl,fl,ff)l,2,3|l)o, 

1^4) = |fl,e,ff)i, 2 ,3|0) OJ 
1 ^ 5 } = ifl,/,s)l,2,3|0) c , 

l^e) = |fl,fl,e)i,2 >3 |0)c, 

I V’t) = |fl,fl,/)i,2,3|0)c. (2) 

Here it is worth noting that in a natural case, the atoms are usually in the same state initially, i.e., the steady state | g). So, 
it is necessary to prepare the initial state \ipi) before implementing the scheme. That is, a population transfer |g)i —> |/)i is 
imperative before the scheme. Fortunately, such operation is not hard to be realized, n pulse, stimulated Raman adiabatic 
passage, large detuning dynamics, and many other techniques are applicable to transfer population from \g)\ to |/)i- Therefore, 
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FIG. 2: Dependence on t/tf of Di/Do and D s /Dq. 






FIG. 4: The time evolution of the populations for states |i/ ) 3 ), I'l/’s), and IV’r) via STIRAP in different cases ((a) and (b) 
are plotted through relation P^ = \(ipm\p\'4’m}\, (c) is plotted through relation P® = |(V , m|4'o)| 2 (m = 1,3,5, 7)): (a) when 
Do = A and tf = 40/A, (b) when Do = A and tf = 80/A, (c) when Do = A and tf = 40/A. 


in subspace V, the interaction Hamiltonian is simplified as (we set A*, = A to be constant coupling coefficients) 

Ho — H a i + H a c , 

H a , = D 1 |V>2)(V'l| + «2|V’4><V’5| + D 3 |V’6><V>7| +H.C., 

H ac = A(|</> 2 > + |V’4) + IV’6»W’3|+tf.C.. (3) 


For the sake of clearness, we will describe two traditional different methods (STIRAP and QZD) to generate W states in brief. 
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FIG. 5: The fidelity of the W state generated via QZD versus the variations of tf. 


A. Based on STIRAP 


For the Hamiltonian in eq. m in the single-excitation subspace, we easily find a dark state 


l*o(*)> 


, Vh > hM _ IV’s) , 

y/Nn Gi G2 G3 A 


(4) 


where Nd = (n ^) 2 + (nj ) 2 + (j ^) 2 + (l ) 2 i s the normalization coefficient. If we choose G 2 ,G 3 ,A Gi at first and make 
sure that the adiabatic condition |(4'o|9t4' n )| |£ n | is satisfied, where |\l/ n ) is the nth eigenstate with nonzero eigenvalue 

the initial state \ipi) = |4'o(0)) will follow I^q^)) closely. Then, we slowly decrease O 2 and G 3 while increase Gi until 
Gi = G 2 = G 3 A at time tf. Accordingly, the dark state becomes |4to(t/)) = -^=(\ipi) + ^ 5 ) + ^ 7 )) which is the W state. 
As shown in Fig. [2j to complete this process, we choose the Rabi frequencies as (we set O 2 = O 3 = O s ) 


Gi — 
G s = 


sinaGo exp[ 


— (t — tp — t f / 2) 2 

t? 


Go exp[ 


— (t + to — tf/ 2) 2 


-l-cos aOo exp[ 


tl 

-(t-t 0 -tf/ 2) 2 

11 


(5) 


where Go is the amplitude and {to,t c } are related parameters. To meet the conditions mentioned above, we choose tana = 1, 
to = 0.15t/, and t c = 0.2 tf. Generally speaking, the adiabatic condition is satisfied better with a relatively large Go because the 
nonzero eigenvalue is proportional to Go. Fig. [3] shows the fidelity of the W state in adiabatic scheme versus the interaction 
time tf and Go. The fidelity F for any target state | ijj) is given through the relation F = |(V’IpI'0)1 5 where p is the density 
operator given through p = i[p, H], As shown in Fig. [3j the fidelity is getting higher with both the increases of Go and tf. It 
seems that when Go x tf > 40, a high-fidelity W state is achievable. That means when Go is large enough, we also can create 
a W state in a short interaction time via adiabatic passage. But further investigation tells us that the system’s evolution is far 
different from adiabatic dark-state evolution with a relatively large Go and a short interaction time tf. Fig. 0(a) shows the 
time-dependent populations for states {|?/>i}, l^s), \tp 7 ), IV^)} when Go = A and tf = 40/A, Fig. 0(b) shows the time-dependent 
populations for states IV’s), \i> 7 ), I 1 / 3 )} when Go = A and tf = 80/A, and Fig. 0(c) shows time-dependent evolution 

of the dark state |4to(t)). The comparison of these three figures draws a result that even with a large Go, a relatively long 
interaction time is still necessary to make sure the controlling parameters change slowly enough to allow adiabatic passage from 
an initial state to a target state. In addition, because a relatively large Go might cause that the RWA is no longer effective for 
the system, and it also means a great population of state \ip 3 ) including a cavity-excited state that makes the system sensitive 
to the cavity photon leakage, it is better to choose a relatively small Go and a long interaction time tf for an adiabatic process. 


B. Based on QZD 

Before we start using QZD to create a three-atom W state, we set G 2 = G 3 = G s and use two orthogonal vectors |/n+) = 
+ |^e)) an< 4 \p~) = — |^e)) to rewrite the Hamiltonian in eq. J3]) as 

H a i = Gi|^ 2 )(^i| + ^=|/x+)((^5:| + (V’rl) 

+ 7§l /x -X(V ,5 l - (VH ) + H.c., 

H ac = A(|V> 2 > + v / 2| M+ »(V> 3 | + H.c.. 


(6) 
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FIG. 6 : The fidelity of the W state in the shortcut scheme versus the interaction time At/ and the detuning A/A in the shortcut 
scheme. 


It is obvious that when the initial state is IV'i), the terms containing |/r-) are negligible because they are decoupled to the time 
evolution of initial state. Then, under the condition <C \/3A (Zeno condition), the subspace V is split into three Zeno 

subspaces according to the degeneracy of eigenvalues of H ac , 


z 0 = (IV’i), |^s), |?M, |<M}, 

Z + = {|0 2 >}, Z- = {\<t> 3 )}. (7) 

where 

I0i> — -^=(-V2|^ 2 > + |m+», 

l<M = -^=(|^2) + Vs\ips) + V2\n+)), 

103) = -^(h/2>-v/3|</3> + V2| M+ )), (8) 

corresponding eigenvalues ei = 0, £2 = a/ 3A, and £3 = — \/3A. Under the Zeno condition, we obtain the effective Hamiltonian 
governing the evolution 


H z = -^|0 1 >(V> 1 | + ^|0 1 >(C| + ^.c., (9) 

where |£) = -^=(| 1 / 5 ) + 1 / 7 )). When Hi and fl s are constant parameters, the general evolution of eq. (9) by solving the 
Schrodinger equation idt |V , (^)) — Hz\4>{t)) in time t is 


10 (*)> 


H 2 + cos fit . , . ifli sin fit.. . 

HW H-7^71—I0U 


+ 


3/f 2 lri/ ' ^ 

cos fit 


3 fi 2 


10 , 


( 10 ) 


where /3 = ^ - n i+ n s ~ g y choosing = (1 ± \/3)Hi and t = tj w the final state becomes | 0 (i/)) = ^ 75 (IV'i) + I'/’s) + IV' 7 )) 
which is the W state. 

In general, the interaction time required in a Zeno scheme is shorter than that in an adiabatic scheme. For example, in 
the present scheme, when we choose relatively large laser pulses, i.e., fli = 0.05A, the interaction time in the Zeno scheme 
is only about tf = 7 t /<5 ~ 35.4/A. However, it is well known that the interaction time should be controlled accurately in a 
scheme via QZD. We plot the fidelity of the W state in Zeno scheme versus the variation in tf in Fig. [5] Here we define 
Sx = x' — x as the deviation of any parameter x, where x' is the actual value and x is the ideal value. It is clear that a 
deviation \Stf/tf\ = 10% causes a reduction about 3% in the fidelity, which shows the scheme is sensitive to the variation of the 
interaction time. In experiment, if we choose a related parameter A = 1GHZ, the required interaction time is tf = 3.54 x 10~ 8 s. 
That means, the experimental researchers should accurately control the interaction time to ensure the deviation in tf is less 
than \5tf\ = 3.54 x 10~ 9 s. That is really a challenge in the current experimental technology. Moreover, known from eq. (1101) . 
the intermediate state \ 4 >i) including atomic-excited states would be greatly populated during the evolution if Hi is too large, 
and that might make the system sensitive to the spontaneous emission. 


III. USING STAP TO FAST GENERATE A W STATE 


Different from the two methods (adiabatic scheme and Zeno scheme) mentioned above, we start from finding a Hamiltonian 
H(t) which is related to Ho{t) to fast generate a W state via STAP. The key point to construct shortcuts for a system governed 











6 


by Ho(t) is to find out a Hamiltonian H(t) which drives the instantaneous eigenstates (I'l'm)} (m = 0,1, 2,3, 4) of Ho(t) exactly. 
Known from Berry’s general transitionless tracking algorithm [Hil ]. the Hamiltonian H(t ) can be reverse engineered from Ho(t). 
And disregarding the effect of phases, the simplest Hamiltonian H(t) is derived in form of 

4 

H(t) = i J2 IftiMOM. (fi) 

m =0 

However, it seems impossible to directly design such Hamiltonian from Ho(t) according to eq. © because the eigenstates 
(|4/ m )} given by solving the intrinsic equation Ho\^ m (t)) = Sm(t)|4' m (t)) are very complex such that mathematically solving 
eq. (HU) seems an outstanding challenge. So here we make a limiting condition -C A to simplify the calculation. 

Under this condition, by substituting the instantaneous eigenstates of Ho in eq. 0 into eq. (1111) . we obtain the 

Hamiltonian that exactly drives the eigenstates of Hz(t) 

4 

H(t) = iY, 

m =0 

= -i0\t){4n\ + H.c., (12) 

where 9 = arctan v ^° 1 , Q a = Q 2 = H 3 , and |f) = -^=(|0s) + | 07 )). Obviously, there is no way to directly apply a pulse between 
the state \ipi) and |£). So, to realize such a Hamiltonian, we need to find out an alternative physically feasible (APF) system 
whose effective Hamiltonian is equivalent to H(t). The model of the APF system is the same as that in Fig. |U (a). The 
difference happens in the atomic level configuration as shown in Fig. [I] (c), each atom also has three levels |/), | g), and |e). 
The transition |/) o |e) is non-resonantly driven by classical field with time-dependent Rabi frequency Cl(t) and detuning A. 
The transition | g) |e) is coupled non-resonantly to the cavity with coupling A and detuning A. Similar to the transformation 
from eq. © to eq. m, the interaction Hamiltonian in the subspace V (the single-excitation subspace for this model is also 
spanned by eq. ©) for the present model can be described as (we set 172 = ^3 = Q 3 and Ai = A 2 = A 3 = A) 

Ho — Hal + Hac + -He, 

Hal = Hl|V’2>(l/’l| + ^=|M+X(0 s| + (0t|) 

+^=Im-X(0s| - (VM) + H.c., 

Haa = Ha C = \{\^2) + V2\^ + )){^\ + H.C., 

H e = A\ip 2 )(ip 2 \ +A\^+)( f j,+ \ +A\n-)(tj,-\. (13) 

The terms including |/r_) are also neglected because they are decoupled to the time evolution of initial state when the initial 
state is set as |0i). Using the eigenstates of H ac to rewrite this Hamiltonian and performing the unitary transformation 
U = e~ l ^ act , we obtain 


l03/jV0i| 


Kt = ni(-^=|^> + -^e i£2t |^> + ^=e 

+fi s (^=l0i) + -±=e i£2t |0 2 ) + ^=e i£3t |</»3»(CI + H.c., 
= A|0l)(0l| + —(102) + |03) ) ((021 + (03|)- 


H, 


(14) 


Consider £ 2 , £3 Hi/\/6, Q s /y/3, we neglect terms containing high oscillating frequencies and terms decoupled to the time 

evolution of initial state. Then we obtain an effective Hamiltonian 

Hz = (-^|0i>(0i| + -^J|0i>(C| +H.c.) 

+A|0i)(0i|. (15) 

1 ) under large detuning condition V ^ 1 , A, we obtain an effective Hamil- 


Then by adiabatically eliminating the state |0i 
tonian 


„ 2|f) 1 | 2 , w | |f7 s | 2 IAWAI 

Heff = -i^H'0i)('0i| — X^HCXCI 

v/ 2 HCH s 


+( 3A IO<V>i| + — 3 ^ 


l0i)<CD- 


(16) 
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When we choose fi s = 0 X and fii = — (here 0 X is a real number), 


H eff = -n\ipi)(ipi\ - n|C)(CI + (-^ICXV’il + h.c.) 

= -n-/ + (-*n|c>(v>i| + h . c .) (17) 

— q2 __ 

where 0 = It is not hard to find, the first term in eq. El) only affects the global phase for the dynamics governed by -ff e //- 
So, in fact, we can directly take off the first term and further simplify the Hamiltonian as H e ff = — + H.c. when we 

pay no attention to the global phase. That means, as long as Cl = 8 and iJ e //(t) = H(t); the Hamiltonian for speeding up the 
adiabatic dark-state evolution governed by Ho under condition \/3A, v^A has been constructed. As we mentioned 

above, to create a three-atom W state by adiabatic dark-state evolution, the Riba frequencies fli and O s can be chosen in the 
form in eq. ©■ Hence, O x is given 


O x = V'SAO 




P 2 


(18) 


If we set l 1 = jj . according to eq. © we can obtain two dimensionless parameters 

t'tf — to — 0.5 tf 
To ’ 

t’tf + to — 0.5 tf 

tc 

Therefore, putting eqs. © and El into eq. El, we obtain 


V i = 
yi = 


(19) 


n 


x 


6\/2AG 2 


tf 


( 20 ) 


where 


G= \ - 


y\Q,iQ ,2 — fliHo(t/ 2 e + cosayie y i) 


20\ + Q.1 


( 21 ) 


is a dimensionless wave function. A brief analysis of G tells us that the amplitude of G is close to 1. That is, the amplitude of 
O x is mainly dominated by v = ■ According to the limited conditions above, we have 


6\/2A 2y/2A 

< V3A t f > 


tf 


6%/2A 


tf 


< V^A t f > 


A 2 

2V2 


( 22 ) 


When A is a constant value, under the premise that the interaction time tf is short, to meet the first condition in eq. (1221) . it 
is better to choose a smaller A, while to meet the second condition, a larger A is required. This is also demonstrated in Fig. 
[ 6 ] which shows the fidelity of the W state in the shortcut scheme versus parameters At/ and A/A. We can find, too small or 
too large A cause a long operation time for the scheme. Then, we choose a set of suitable parameters {A = 3A ,tf = 35/A} for 
the scheme. As shown in Fig. [7] (a), with these parameters, we can achieve a perfect populations transfer after very slightly 
correcting a related parameter (Hi —> 1.04Hi) by numerical simulation. The parameter should be slightly corrected because 
speeding up the evolution needs to slightly broke the Zeno condition, and that cause a slight failure of the approximation. We 
plot the time evolution of the populations in intermediate states l'i/ 3 ) and |</i) in Fig. |7j(b) to prove this operation is necessary. 
In the figure, the state \ 1 jj 3 } should have been neglected is slightly populated (the Zeno condition is slightly broken) while the 
state |0i) is negligible (the second condition in eq. (1221) is fulfilled). We give a comparison of the fidelities via these three 
different methods in Fig. [7](c). Contrasting with the Zeno method, the advantage of the present shortcut method is obvious: 
the shortcut scheme is more robust against operational imperfection than the Zeno one, especially, it is not necessary to control 
the interaction time accurately. As demonstrated in Fig. [H] the fidelity almost keeps unchanging with the variation ST, where 
T = 40/A is the total operation time chosen to complete the scheme, and a deviation \ 8 v/v\ = 5% which means the variation 
in the amplitude of O x only causes a reduction about 1 % in the fidelity. 

Now, we will check the robustness of the shortcut scheme against possible mechanisms of decoherence. The evolution of the 
system can be modeled by a master equation in Lindblad form when the decoherence is considered, 

p = i[p, H 0 ] + J2^ Lk P L l - \( L l L kP + pL{L k )\, 

k 


(23) 
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FIG. 7: (a) Time evolution of the populations for the states IV’i): IV’s) an< A l^/r) with {tf = 35/A, 5 = 3A}. (b) Time evolution 
of the populations for the intermediate states \tp 3 ) and |</>i) with {tf = 35/A, S = 3A}. (c) The comparison between the fidelities 
of the three schemes, the blue solid curve representing the shortcut scheme is plotted with {tf = 35/A, S = 3A}, the red dash-dot 
curve representing the adiabatic scheme is plotted with {tf = 80/A, fio = A}, and the green dash curve representing the Zeno 
scheme is plotted with {tf ~ 35.4/A, fli = 0.05A}. 


where p is the density operator for the whole system and L k are the Lindblad operators. For the shortcut scheme, there are 
seven lindblad operators governing the dissipation: 

Li = V7T|/)i(e|, L% = •v/ 72 |s i )i(e|, 

L 3 = V73|/>2<e|, Ll = 

Ll = v / 7sl/)3(e|, Lg = ^\g) 3 {e\, 

L 7 = V^a, (24) 

where 7 m (m = 1, 2, • • • , 6 ) are the atomic spontaneous emissions and re is the cavity decay. We set 7 m = 7/2 for simplicity. 
Then by numerically solving the master equation in eq. (1241) . we plot the fidelity of the W state in the shortcut scheme versus 
7 /A and re/A in Fig. [21(a). We can find that the shortcut scheme is more sensitive to the cavity decay than atomic spontaneous 
emissions with parameters {A = 3A, tf = 35/A}. The reason has been mentioned above that with this set of parameters, the 
Zeno condition is not satisfied faultlessly. So the states |</> 2 ) and \(j} 3 ) containing cavity-excited state l'i/ 3 ) are populated in a 
certain extent during the evolution. Furthermore, the parameters can be selected properly to restrain the cavity decay in the 
experiment according to eq. (12211 . For example, when we choose {A = A , tf = 35}, the Zeno condition can be satisfied well. 
We plot Fig. [2] (b) depicting the fidelity of the W state governed by the APF Hamiltonian versus re/A and 7 /A when (A = A, 
tf = 35}. It shows that the influence of cavity decay is restrained with these parameters. However, it is without doubt that 
the scheme is robust because the fidelity decreases slowly and even when 7 = re = 0.1 A, we still can create a W state with a 
high fidelity 91.12%. 

This scheme can be easily generalized to generate IV-atom W states. We assume N A-type atoms are trapped in a cavity. 
For the original Hamiltonian, the atomic level configuration of each atom is the same as that in Fig. [l](a), and for the APF 
Hamiltonian, the atomic level configuration of each atom is the same as that in Fig. [T] (c). Suppose that the N atoms ‘see’ 
the same field, and spatial separation size of these atoms is much bigger than the wavelength of the emitted radiation, so the 
atomic dipole-dipole interaction can be omitted. In this case, the interaction Hamiltonian for the original Hamiltonian reads 

N 

H? = ]T fi fc (t)|e)fc</| + \ k a\e)k(g\ + H.c., (25) 

k= 1 


and the APF Hamiltonian reads 


N 

H? = ^2 A l e M e l + ^fc(*)|e>fc</| + X k a\e)k{g\ + H.c.. 

k= 1 


(26) 
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FIG. 8: The fidelity of the W state via STAP versus the variations of T and v. 


We consider that the initial state of the system is in \ f, g,g, ■ ■ ■ ,g) 1 , 2 , 3 ,--- ,jv| 0) c . For the atom 1, the classical field drives the 

transition resonantly between the level |/)i and |e)i with the Rabi frequency fii. Then, atom 1 will emit a photon which 

will be absorbed by one of the other AT — 1 atoms with the same probability when Ai = A 2 = ■ • • = Xn = A. Therefore, 
the excited process of N — 1 atoms (in this part, the “N — 1 atoms” means the atoms except the atom 1) can be described 

by the state |$ e ) = -j={\e,g,g,-■ ■ , g) + \g,e,g, - ,g) + \g,g,e,--- , g) H-b | g,g,g,--- , e» 2 , 3 , 4 ,- ,n- Then, by setting 

O 2 = O 3 = ■ ■ ■ = fijv = O s , the classical fields will drive the state |'h e ) to (|/, g,g, ■ ■ ■ ,g) + |g,/,<?,••• ,g) + 

\g, g, /, • ■ • ,g) + •■■ + \g,g,g,-" , /)) 2 , 3 , 4 ,-. ,n- Hence, the single-excitation subspace could be spanned by 

IVh) = \f,g,g,-‘- ,3>i,2,3,-.,jv|0)c 

IX 2 ) = je,<?,<?,••■ , 5 )i, 2,3,-.. ,jv|0) c , 

IV’3> = \g,g,g, ■■■ ,s)i,2 1 3,- 1 jv|i) c , 

lM> = |5>l|®e>|0)c, 

10 = lfl)i|®/>|0>c. (27) 

Meanwhile, the Hamiltonian in the single-excitation subspace can be written as 

H f = n^XlM + n s \g}((\ + A(|^ 2 ) + 0]V^T|m»<V>3| + H.C.. (28) 

Similarly, the APF Hamiltonian in the single-excitation subspace is 

H? = Ad^X^I + l^eX^el) 

+[f 2 i|V> 2 )(V-i| + n.|/*)(C| + A(|^ 2 > + T|a*»<^I + H ^}- (29) 

Obviously, the Hamiltonians in eqs. (EU and (O are in the same form with those in eqs. © and respectively. Therefore, 
similar as above, under the condition £7i, fi s , fX, f2 s ^ |e±| and fli, <C A, where e± = iy/N A are the nonzero 

eigenvalues of H^ c = X(\tp 2 + y/N — l\g))(tp 3 \ + H.c., and H^ will be approximated as Hz = — + 

7 ^l<^i>(CI +H.c. and Hz = _ (Jv ^| ni|2 1 (ipi h 10(Cl + ( lQs 10(Xi 1 +H.c.), respectively. By setting fX = ~^== 
and 0.2 = ^3 = ■ • ■ = On = O x , the effective Hamiltonian which is equivalent to the counter-diabatic driving Hamiltonian of 
the original Hamiltonian will be achieved. Then, the shortcut can be constructed and the A-qubit W states can be rapidly 
generated. 

In a real experiment, the cesium atoms which have been cooled and trapped in a small optical cavity in the strong-coupling 
regime [60l.[61j can be used in this scheme. On the other hand, a set of cavity QED parameters (A, 7 , k)/2-j: = (750, 2.62, 3.5)MHz 
is predicted to be available in an optical cavity [62|. With these parameters, the fidelity of the W state in the shortcut scheme 
is 99.01%. 


IV. CONCLUSION 

In this paper, we have proposed a scheme to fast generate W states via transitionless-based shortcuts. In order to highlight 
the advantages of the present scheme, we have described two similar schemes based on STIRAP and QZD. The comparison 
among these three schemes demonstrates that the shortcut scheme is faster than the adiabatic one, and more robust against 
operational imperfection than the Zeno one. Numerical investigation also demonstrates the present scheme is robust against 
the decoherence caused by both atomic spontaneous emission and photon leakage. When it comes to the generation of Af-atom 
W states, the only change is setting fX = — -f=== and ^2 = fl 3 = • • • = = fire. Known from ref. [63], the Hamiltonian for a 

system with three four-level atoms trapped in a cavity also can be approximated into an effective Hamiltonian in form of eq. ©■ 
For a similar model described in ref. [6J], if the atomic transitions are non-resonant, one can also obtain an effective Hamiltonian 
in form of eq. m- That means, with the same method in section III, a multi-qubit singlet state also can be fast generated. In 
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FIG. 9: Dependences on k/A and 7 /A of the fidelity of the W state governed by the APF Hamiltonian when (a) tf = 35/A and 
A = 3A; (b) t f = 35/A and A = A. 


addition, the shortcut method might also show its glamour in other fields, for example, fast transfer of entanglement [32, 65j] 
That demonstrates the present method has a wide rang of application in quantum information processing. This might lead 
to a useful step toward realizing fast and noise-resistant quantum information processing for multi-qubit systems in current 
technology. 
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